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Abstract 

We study the evolution of tensor metric fluctuations in a class of non-singular mod- 
els based on the string effective action, by including in the perturbation equation the 
higher-derivative and loop corrections needed to regularise the background solutions. 
We discuss the effects of such higher-order corrections on the final graviton spectrum, 
and we compare the results of analytical and numerical computations. 
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1 Introduction 



Recent studies, including strong coupling corrections to the string effective action, have pro- 
vided the pre-Big Bang scenario |[[ |2| with a number of promising models for the transition 
from the growing to the decreasing curvature regime || [|, |], H . However, a crucial question 
remains: is it possible to test this cosmological scenario and, in particular, to distinguish it 
from other, more conventional descriptions of the very early Universe? 

The answer to this essential question is in principle positive, since the pre-Big Bang scenario 
incorporates a dilaton field which modifies the standard inflationary kinematics and which, 
through its non-minimal coupling, also directly modifies the evolution of perturbations (of the 
metric, and of the other background fields) 0. It is well known, on the other hand, that the 
transition from an inflationary period to the usual Friedman-Robertson- Walker (FRW) phase 
amplifies tensor metric perturbations, and is associated with the production of a stochastic 
background of relic gravitational waves || [|. 

Such a primordial background decouples from matter immediately below the Planck scale, 
unlike the electromagnetic radiation which underwent a complicated history until recombi- 
nation, and has been transmitted almost unperturbed down to our present epoch. As a 
consequence, its present spectrum should be a faithful portrait of the very early universe [10|, 



thus opening a window for the observation of processes occurring near to the Planck scale, 
and for discriminating among various models of primordial evolution. 

In the context of string cosmology, several mechanisms may contribute to the generation 
of gravitational waves from the initial vacuum state. One contribution comes from the process 
of dynamical dimensional reduction [IT, [T2], [TJ|, when the internal dimensions shrink down 



to a final compactification scale. Another contribution is due to the time-dependence of the 
dilaton field - and thus of the effective gravitational coupling constant - during the pre-Big 
Bang phase [Tj|, ITS]. Finally, there is the usual contribution arising from the accelerated 



expansion of the external three-dimensional space. Both the dilaton and the full, higher- 
dimensional metric thus contribute to the external "pump" field which is responsible for 
the parametric amplification of tensor metric fluctuations, normalised to an initial vacuum 
fluctuation spectrum. 

With respect to the standard inflationary scenario, the amplification of tensor perturba- 
tions in the pre-Big Bang scenario is strongly enhanced for large comoving wavenumber k 
TEL |T4}1, in such a way that the produced background of gravitational waves could in principle 



be detected in the near future by various planned experiments [EJ, [lTj, I7|. Arguments in favour 
of this conclusion arise from the particular kinematics of string cosmology. First, since the 
Hubble parameter is increasing in the pre-Big Bang phase, higher frequency modes will cross 
the horizon at higher values of \H\, which implies that the amplitude of the spectrum increases 
with frequency. In this context, the comoving amplitude of perturbations may grow even out- 
side the horizon |TB|, [TTj, ^TJ , instead of being frozen as happens in standard inflation. Second, 
the peak amplitude of the spectrum (naturally located in the high frequency regime, near the 
end point of the spectrum) is normalised so as to match the string scale ^TJ, i.e. nearly 
eight orders of magnitude above the high frequency normalisation of the standard scenario 



22]. Hence, the pre-Big Bang phase could lead to a very efficient production of gravitons in 



the sensitivity range of present gravitational antenna. In fact such a relic background could 
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be in principle be detected by the second (planned) generation of interferometric detectors 



23 



To date, extended studies of tensor perturbations have been mostly performed in the con- 
text of the tree-level effective action, with and without the possible contribution of a high 
curvature stringy phase. In all cases the background cosmologies that have been investigated 
have been singular, i.e. the curvature and dilaton have become singular at some fixed proper 
time. There have been a number of suggested ways to tackle the singularity problem, con- 
sidering also anisotropic |25|, [26|] and inhomogeneous |27, 28] backgrounds, or the presence 
of a non-local dilaton potential 0. However, if we confine our attention to homogeneous 
and isotropic metrics, and to a local potential, one of the most promising approaches to the 
graceful exit problem |29j suggests that the curvature singularities may be cured by adding 
higher-order corrections to the string effective action (see for instance [31], ^2[ ||, £|, ||). 



These include both higher-derivative a' corrections, due to finite-size effects, and quantum 
corrections, due to the loop expansion in the string coupling parameter. In the context of 
cosmological perturbation theory, it is known that such corrections may induce an addi- 
tional amplification when the curvature becomes large in string units. The aim of this paper, 
therefore, is to discuss the evolution and the spectrum of tensor perturbations in the context of 
a class of non-singular cosmological backgrounds, by taking into account (in the perturbation 
equation) the full contribution of those corrections that are responsible for the regularisation 
of the background evolution. In this way it will be possible (for the first time, to the best of our 
knowledge) to follow the complete evolution of perturbations, from the initial vacuum down 
to the present time, through a continuous numerical integration of the linearised perturbation 
equations. 

The paper is organised as follows. In Section 2 we recall the general form of the first order 
ol corrections arising in the context of the massless bosonic sector of the low-energy heterotic 
string action. We extend them by including possible one- and two-loop quantum corrections, 
and we write down the corresponding background equations. In Section 3 we compute the 
linearised tensor perturbation equation, including such higher-order corrections. In Section 4 
we present a general discussion of the evolution of perturbations in the presence of higher-order 
curvature corrections. We then check the main points of our discussion through a numerical 
integration of the perturbation equations for regular backgrounds, by using an appropriate 
dilaton potential, or particle production effects, to stabilise the dilaton in the post-Big Bang 
era. We find that the general effect of the higher-order corrections is to flatten the slope of 
the high-frequency part of the spectrum, but the effect is probably too small - in the class of 
models discussed in this paper - to make the graviton background "visible" to the planned 
advanced detectors. The main results of this paper are finally summarised in Section 5. 



2 Effective action and background equations 



In the context of the pre-Big Bang scenario, our present FRW Universe is expected to emerge 
from the string perturbative vacuum as a consequence of some process - for instance, plane 
wave collisions |34j - triggering the gravitational collapse of a sufficiently large portion of 
spacetime [35], The initial evolution is driven by the kinetic energy of the dilaton, and can be 
described, in the string frame, as inflationary expansion with growing coupling and curvature. 
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After reaching a maximal scale controlled by the fundamental string length As, \H\ ~ O^X^ 1 ), 
it is then expected that the Universe is smoothly connected to the FRW regime, with a constant 
dilaton field. 

In the early stage of the dilaton-driven inflationary (DDI) period, strings are propagating 
in a background of small curvature, and the fields are weakly coupled, so that the cosmolog- 
ical evolution is consistently described by the lowest order string effective action which, by 
assuming static internal dimensions, can be written as: 



r = 1 J d\^~g {e"* [R + (d,4>) 2 ] + C c ) . (1) 



Here a' = A|, we have adopted the convention (— , +, +, +), -R^ Ap = d p T^ x + . . ., R up = R x uXp 
and set our units such that h = c = 167rG = 1. Following ||, we have introduced an additional 
lagrangian C c allowing for the inclusion of higher-order corrections to the tree-level effective 
action, a non-perturbative potential - yet to be determined - and/or the backreaction of the 
produced radiation. Considering this additional lagrangian as an independent source we shall 
assume that the associated energy-momentum tensor is diagonal, homogeneous and isotropic, 
and can be written in the perfect fluid form as T£ = diag(— p c ,p c ,p c ,p c ). 

The resulting equations of motion for the background are obtained from the variation of 
Eq. (H) with respect to the metric and the dilaton field: 

QH 2 — QH<j) + <j) 2 = eV, (2) 

4<j)H - AH - QH 2 - <p 2 + 20 = e p c , (3) 

QH + 12H 2 + 2 - 20 - 6<j)H = e A^, (4) 

p c + 3H( Pc + Pc ) = -A^0. (5) 

A dot represents the differentiation with respect to proper time t and H = a/a is the Hubble 
parameter. arises from the variation of the additional lagrangian with respect to <fi. Finally, 
Eq. (|) is the conservation equation, a direct consequence of Eq. (^)-Eq. (|). 

To lowest order, the homogeneous and isotropic version of the pre-Big Bang scenario, 
without non-local dilaton potential, inevitably faces a curvature and dilaton singularity at 
the end of the DDI phase. However, a singular behaviour is often a late manifestation of the 
breakdown of the description, extended beyond its domain of validity. Hence, the low-energy 
dynamics should be supplemented by corrections, in order to give a reliable description of 
the high curvature regime, expected to be crossed before reaching the FRW universe. Such 
corrections are generally of two types: tree-level a 1 corrections, referring to the finite size of 
the string, and quantum corrections, resulting from a more conventional loop expansion in 
the string coupling parameter. 

It has been shown in [|32| that the higher curvature a' terms, arising from the expansion 
around the point-particle limit, may eventually stabilise the evolution of the Universe into 
a de-Sitter like regime of constant curvature H ~ 0{X^ 1 ) and linearly growing (in cosmic 
time) dilaton. In this paper we will consider the most general tree-level a' corrections to the 
gravi-dilaton string action up to fourth-order in derivatives (see for a detailed analysis), 



which can be written in the form 36 : 



C a > = a'Xoe^ { Cl R 2 GB + c 2 G"^<R0 + c 3 D0(^0) 2 + c 4 («9 M 0) 4 } , (6) 
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ci = -1, c 2 + 2(c 3 + c 4 ) = -2ci. (7) 

The constraint Eq. (|7|) on the coefficients is required so that the action reproduces the usual 
string scattering amplitudes |5B|. The parameter X Q allows us to move between different 



string models, and we will later use A = — 1/4 to agree with previous studies of the heterotic 
string Finally, R 2 GB = R^xpR^ ~ ^R^R^ + R 2 is the Gauss-Bonnet combination, 

guaranteeing the absence of derivatives higher than two in the equations of motion. 

Another conventional expansion of the string effective action is controlled by the string 
coupling parameter, g 2 tring = e®. Initially, g 2 string <C 1. However, the growth of the dilaton 
during the DDI era implies that the late evolution of the background could be dominated by 
quantum loop corrections, which will eventually force the universe to escape from the fixed 
point determined by the a' corrections fl32| |, for a smooth connection to the post-Big Bang 



regime. Unfortunately, there is as yet no definitive calculation of the full loop expansion in 
string theory, so we are left to speculate on plausible terms that will eventually make up the 
loop contribution. Multiplying each term of the tree-level a' correction by a suitable power 
of the string coupling is the approach which has already met with some success p, ||, 
and which we shall adopt also in this paper. Since the quantum corrections are not formally 
derived from an explicit computation, we shall allow for different coefficients di (e^) at one-loop 
(two-loop) in the string coupling, replacing the coefficients q of the tree-level action, where 
di and are not necessarily subject to the constraint Eq. (^). With such assumptions, the 
effective lagrangian density C c of Eq. (H), which is in general the sum of the tree- level a' and 
loop corrections, C c = C a i + C g , in the present case will take the form 

C c = CM + Ae^Caidi) + Be 2 *£ a ,( ei ), (8) 

where C a '(ci) is given by Eq. (||), with the constant parameters A and B actually controlling 
the onset of the loop corrections. 

The very similarity between the first order a' and loop corrections enables a simple form 
for their contributions to the background equations of motion. From now on, we introduce 
an additional parameter n in the exponential of the correction C c , such that a' corrections 
correspond to n — 0, whereas one- or two- loop corrections correspond to n — 1 and n = 2 
respectively. For instance, using the notation {. • .}i to represent the one- loop contribution 
(n = 1 inside the brackets), we easily get the source terms for the system Eq. (Q)-Eq. (f|): 

Pc = {Pc}o + A{ Pc } 1 + B{p c } 2 , (9) 
Pc = {p c }o + A{p c } 1 + B{ Pc } 2 , (10) 
= {A^o + ^VK + ^A^h, (11) 

and the terms in brackets are given respectively by: 

{p c } n = a / A o 0e (n - 1 ^{-24c 1 (n - 1)H 3 + 9c 2 <pH 2 + Gc 3 <j) 2 H - [c 3 (n - 1) - 3c 4 ] 3 }, 



{p c } n = a 'X e {n ~ 1 ^{ [8ci(n - l) 2 - 3c 2 ] <p 2 H 2 - 2c 2 (n - l)^H - 2o 



; 30 2 <; 



+8d(n - l)(f)H 2 + 16ci(n - l)<f>H 
-4c2H>H - 2c 2 2 ij + [c 4 - c 3 (n - 1)] 4 }, 



{Mr 



a 



'X e {n ^{2Ac 1 (n-l)H 



in 



l)> d + 4(n- 1 

3 



3c 4 2 [ty + 4<pH +(n-l)<j) 2 
(H + 3H 2 ) - 120H]} 

#]}. 



-3c 2 [24>H 2 + 6<pH 3 + 4<pHH + {n - V : 2 

Here, of course, for n = 1 the coefficients q are replaced by c?j, while for n = 2 the c$ are 
replaced by e«. Note that for n = 0, Aq = —1/4, C4 = — ci = 1, one easily recovers the a' 



corrections used in [32, 33 



In general, the combination of tree-level a' and quantum loop corrections does not lead 
to a fixed value for the dilaton in a finite amount of time. Although a non-perturbative 
(supersymmetry breaking) potential [[37] is expected to stabilise the dilaton, here we shall 
follow p| by assuming the dilaton is frozen out by radiation, after the transition to the post- 
Big Bang regime. To this aim we will introduce "by hand" the presence of radiation, coupled 
to the dilaton and satisfying the conservation equation 



p + AHp- h 



0. 



(12) 



where represents the decay width of the dilaton. 

The solutions of the system of equations Eq. (0) - Eq. (||), including the radiation and the 
quantum corrections, provide the cosmological gravi-dilaton background in which we will study 
the propagation of tensor perturbations. We will consider a spatially flat FRW manifold, and 
we will perturb the above equations keeping the dilaton and all sources fixed, 5(f) = 0, 5T^ = 0. 
To lowest order, the tensor fluctuations will obey the usual string frame perturbation equation 



14], including the non- minimal coupling to a time-dependent dilaton. On the other hand, 



the inclusion of higher-order corrections, needed to regularise the background, will inevitably 
lead to a modification of such a perturbation equation (see P|L |39], HO, Hl[ for initial studies 



taking into account higher curvature contributions to the evolution of perturbations) 



3 Tensor perturbations 

It is well known that gravitational waves, arising from linearised tensor perturbations, do not 
couple to pressure and energy density, and hence do not contribute to classical gravitational 
instabilities. Nevertheless, they are of interest as a specific signature of pre-Big Bang cos- 
mologies since, as stressed in a number of papers [0, [16|, |14||, the production of high-frequency 
gravitons is strongly enhanced, in string cosmology, with respect to the standard inflationary 
scenario. In this study we will restrict our attention to tensor metric perturbations around a 
d = 3 spatially flat background, parameterised by the transverse, trace-free variable h^ u , 

g^^g^ + h^ with v u h^ = = h£, (13) 

where V u denotes covariant differentiation with respect to the background metric. Note also 
that the indices of hu V are raised or lowered with the unperturbed metric, h v = g up h fip . 
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3.1 Linearised equations 



The linear evolution equation for the metric fluctuations can be obtained by perturbing the 
action Eq. ([!]) up to terms quadratic in the tensor variable h^ v : 



+A J d A x [{5 {2) £} + A{5®C}i + B{5 {2) £} 2 
where the integrand of the n-correction is given by: 



(14) 



>-i) 



^ 2 \ Cl ^-gR 



6 (2) 



C2V-9 



RK» _ -g^R 



+5^ 



C3V-9 1 



□0(<9 M 0) 2 ) + <f< 2 > (ctyf^id^) 



(15) 



Following fl33fl , it is convenient to use the synchronous gauge where 
fi'oo = — 1) 9oi — 0, Qij = a 2 Sij, 



h 



00 



0. 



h 



0. 



g l3 h 



and djh^ 



0. 



(16) 



This enables us to write the perturbed action as a quadratic form depending on the first 
and second derivatives of the symmetric, trace- free matrix h = h % -, with time-dependent 
coefficients fixed by the background fields a(t), <p(t). Introducing the matrix notation Tr h 2 = 
h-'hj, and using V 2 = S^didj for the Laplace operator in a flat space, we have 



<j(2) r = - [d 4 x a 3 e^ Tr 
a' J 



Itf 3H 2 
V 2 



h 2 



• S ■ 1 V 2 ' 
-hh - AHhh - -h 2 + -h--h 
4 4 a 2 



+A / d 4 x Tr [{®} + A{0} 1 + B{®} 2 ] 
where the integrand of the n-correction is given by: 



(17) 



{®}n 



a 3 e (»-U*< 



-6 Cl H 2 (H + H 2 ) - l -j> 2 (3c 2 H 2 + c 3 4> + 3c 3 0H + 



ciH + 7ci# 2 + -on, 



h 2 - A Cl H 2 hh - 2 Cl Hhh 



%c x H{H + 2H 2 ) + c 2 <j) 2 H + -c 3 3 



. V 2 • 
hh + 2c\h—-h 



+ 



ciH 2 + ciH + -c 2 2 



V 2 -V 2 -V 2 ' 

h—h + AciHh—h + 2c x h—h 



Here, and in what follows, the coefficients q are to be replaced by d{ and for n = 1 and 
n = 2, respectively. Again, for n = 0, X = —1/4, C4 = —c\ = 1, one recovers the first order 
a' corrections discussed in [33|. 



We can now integrate by parts all the terms with more than two partial derivatives acting 
on h, as well as the terms in hh and hh. This will drop all terms with more than two derivatives 
thanks to the Gauss-Bonnet combination, leading to 



£< 2 >r 



a' J 



(fx a 3 e-* Tr 



+-h 2 + -h—h 
4 4 a 2 



+A J d 4 x Tr [{®} + A{®} 1 + B{®} 2 ] 



where we have 

{®}n 



aV n -^|-^ 2 ([8 Cl (n-l) 2 -3 C2 



Z.2 tt2 



H z - 2co(n - 1)6 H 



+ [c 4 - c 3 (n - 1)] 4 - 2c 3 2 + 8ci(n - 1)6H 2 
+16ci(n - l)j>H [H + H 2 ] - 4c 2 00if - 2c 2 2 #) 



V 2 
+h—h 

a z 



4> 2 + d(n- 1)0 



H + H 2 

ci(n - l) 2 + \c 2 



+h 2 ci(n- l)(f)H - -c 2 fr 



We stress that all tree-level a' and quantum loop corrections disappear in the limit of a con- 
stant dilaton field, since in that case the contribution of the Gauss-Bonnet term corresponds 
a topological invariant. This ensures that the resulting perturbation equation reduces to the 
Klein-Gordon form, typical of a minimally coupled scalar, as we enter the frozen dilaton, 
post-Big Bang regime. 

The rather complicated form of Eq. flT8| ) can be simplified as the coefficient of the h 2 term 
vanishes identically, thanks to the background equation Eq. (|3]). By decomposing the matrix 
h l j into the two physical polarisation modes of tensor perturbations, h + and h x , 



Trh 2 = h)h{ = 2(h 2 + + hl), 
we can finally write the action, for each polarisation mode h(x,t), as 

5 {2) T h = i-fAaV^ 2 (l + {®} + A{©}i + B{ffi} s 
2a J L 



(19) 



with 



{©}n 
{®}n 



c/A e n * 



+h—h (1 + {®} + A{®} 1 + B{®} 2 ) , 



4 Cl (n-l)0#-^c 2 2 
4ci(n - 1) 2 2 + \c 2 4> 2 + 4ci(n - 1) 
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(20) 



where h is now a scalar variable standing for either one of the two polarisation amplitudes 
h + , h x . The variation of the action with respect to h gives then the modified perturbation 
equation: 

= [1 + {©} + A{©}! + 5{©} 2 ] /i 

+ [3H - (j) + {®} + A{®} 1 + B{®}^ 

- [i + {e} + A{eh + B{e} 2 ] ^-h, 



(21) 



with 



{©}, 

{®}r 



a'A e n * \ 4ci(n - 1)0# - -c 2 <^ 



a'A e n< ^ 4c x (n - 1) (<j>H + - c 2 



-c 2 (yi — 1, 



a'Aoe 



n<f> 



+ 12 Cl (n - l)0iT + ( 4ci(n - l) 2 - -c 2 



4 Cl (n - 1) 2 2 + -c 2 2 + 4ci(n - 1) 



2 H 



Eq. (pi]) controls the time evolution of the Fourier components of the two polarisation 
modes, and is valid even during the high curvature regime connecting the pre-Big Bang branch 
to our present FRW universe, since it incorporates the contributions of the higher-order cor- 
rections needed to regularise the background evolution. 



3.2 Quantisation in conformal time 



We shall now focus on the generation of gravitational waves in the context of the pre-Big Bang 
scenario. To normalise the graviton spectrum to the quantum fluctuations of the vacuum we 
need the canonical variable that diagonalises the perturbed action, and that represents in this 
case the normal modes of tensor oscillations of our gravi-dilaton background. To this purpose 
we note that introducing the conformal time coordinate rj, defined by a = dt/dr], the action 
Eq. (p0|) can be re-written 



S®r h = / d 3 xd V {z 2 ti 2 + y 2 hV 2 h\ , 

where a prime denotes differentiation with respect to the conformal time, and 

y 2 ^) = e -^[a 2 + {©} + A{©} 1 + J B{®} 2 ] ) 
z 1 (r 1 ) = e-* fa 2 + {®} + A{®} x + B{®} 2 ] , 



with 



{©}r 
{®}r 



a 



'V* 4 Cl l 



n 



a'A e n ^4ci(7i- 1 



a'(j>' 



+ lc 2 /2 + 4 Cl (n-l) U" 



;c 2 ( 



a'(f)' s 



(22) 

(23) 
(24) 



Introducing the canonical variable ip = zh enables us to diagonalise the kinetic part of the 
action: 

S®T H = ±J d 3 xd v jv/ 2 + y^ 2 + f^Vv} • (25) 

For each Fourier mode, ipk = zhk, we can now obtain from the perturbed action a linearised 
wave equation in terms of the eigenstates of the Laplace operator, V 2 ipk = —k 2 ipk which takes 
the explicit form: 

V£ + {^-Vfo)}v* = 0, V( V ) = ^. (26) 

This linear perturbation equation is the major result of the first part of this paper, since it 
encodes (through Eq. (p3|) - Eq. (|24]) ) the full contribution of those higher-order corrections 
which can be used to regularise the background evolution. We note that the above equation 
defines an effective time-dependent mode, ky/z, whereas the effective potential is determined 
by the evolution of the generalised pump field z{rf). These new features of our generalised 
perturbation equation will be discussed in detail in the next section. 

The graviton spectrum of pre-Big Bang models generated in the context of the lowest-order 
perturbation equation, with and without the contribution of a high curvature string phase, 
has been discussed extensively in the literature (see for example @, |2(| [17], [42], |43|, [44]]). In the 



remaining part of this section we shall briefly recall some general results of these studies, so 
as to highlight the modifications induced by the higher-order corrections in the perturbation 
equation, and in the spectrum of the relic pre-Big Bang gravitons. 

In the extremely weak coupling and low curvature regime, the tree-level solutions of the 
pre-Big Bang scenario are fully adequate to describe, asymptotically, the initial evolution of 
the quantum fluctuations of the gravi-dilaton background. In normalising the fluctuations in 
the asymptotic past (r] — > — oo), we shall thus neglect any correction to the low-energy action, 
considering the standard perturbation equation for the Fourier mode i[>k(v)'- 



€ + k 2 



r(r-l) 
(—ki]) 2 



ij k = 0. (27) 



Such an equation describes the evolution of a classical harmonic oscillator undergoing a para- 
metric amplification driven by the effective potential V = w hh £ ~ (— 7]) T . By recalling 
the (conformal time) expression of the tree-level solutions 0, 

a(v) = (-VT, ot = (1 - V3)/2, (28) 

<f ) ( v ) = - 1 ln(- V ), 7=V3, (29) 

we obtain the external pump field £ = ae~^ 2 = (-r/) T , r = 1/2. Hence, the power-law type 
of the effective potential leads to two distinct regimes for the evolution equation Eq. ([27]). 

Sub-horizon modes (\krj\ ^> 1) are freely oscillating, and no quantum particle production 
occurs. However, the non-adiabatic behaviour of modes having left the Hubble radius (\kr]\ 
1) in the dilaton-driven phase, are described asymptotically by the general solution 



fe(,) = * + a.jT ( ^, (30) 
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where r\ exit = A; -1 , and h k = tp k £ _1 is the comoving amplitude of tensor perturbation for 
the given mode (A k and B k indicate the coefficients of the growing and decaying solutions, 
respectively). The amplitude of super- horizon modes is thus frozen, modulo logarithmic cor- 
rections. The overlapping of these two extreme behaviours is given by the exact solution to 



the Bessel-type equation Eq. (27), which can be expressed in a linear combination of the first- 



and second-kind Hankel functions: 

Mv) = (-ki]) 1/2 { Kl Hjp{-kri) + K 2 Hj?\-kri)} , (31) 

where v — |1 — 2r|/2 = encompasses the background dynamics. 

In the asymptotic past (i] — > — oo) any given mode k is well inside the Hubble radius (as 
H^ 1 blows up as we go back in time in pre-Big Bang models), and the two types of Hankel 
functions H^ 1 ' 2 ' (k\r]\) are oscillating, corresponding to negative and positive frequency modes 
respectively. It is therefore possible to normalise tensor perturbations to a spectrum of quan- 
tum fluctuations of the initial background. Indeed, the annihilation and creation operators 
resulting from the expansion of the h^u field, with ip{k,rf) = a k ty + a\,^/*, obey canonical 
commutation relations [a k , a k >] = [a' k , a\,] = and [a k , a k ,} = 5 kk >. For the normalisation to an 
initial vacuum state, we choose to restrict to positive energy states, ip k — > e~ lkri , hence 

Kl = and k 2 = y / fe-^( 1+2 ^. 

The effective potential of the perturbation equation is expected to grow monotonically 
during the initial dilaton-driven regime, up to a maximal value reached during the string 
phase, while it is expected to decrease rapidly to zero at the onset of the standard, radiation 
dominated phase. Emerging from the asymptotic past, a typical mode k will thus first oscillate 
inside the horizon and then progressively feel the potential barrier. The process of quantum 
particle production will take place and carry on until the background enters the standard 
FRW- regime. After that, the effective potential identically vanishes, and Eq. fl2~T|) yields the 
free oscillating solution: 

i, k = -L {c+(k)e- ik * + c_(k)e ik ^ V - +oo, (32) 



2k 

with incoming and outgoing waves. The mean number of produced gravitons is obtained from 
the square modulus of the Bogoliubov coefficient, \c_(k)\ 2 , using the matching of the tensor 
perturbation variable h and h' at the onset of the FRW phase. 

The low-energy perturbation equation may allow for estimating a lower bound on the spec- 



trum of the produced gravitons [ 46[ . For the tree-level solution with r = 1/2, by defining the 



spectral energy density per logarithmic interval of frequency, p gw (k), one obtains in particular 

0. 

where k\ = refers to the maximal height of the effective potential barrier, l^(^i)! 1 ^ 2 — 

l^il -1 - We recall that the cubic slope of the spectrum is a direct consequence of the evolution 
of the background during the DDI branch. In [ 17 , pE3| ], the authors argued that the 



contribution of a high curvature string phase would lead in general to a reduction of the slope 
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of the spectrum of relic gravitational waves. This was explicitly confirmed for a model based 
only on the tree- level a' corrections [|33|. However, we stress that the investigation of the 



evolution of high frequency modes, hitting the potential barrier during the high curvature 
regime, needs to be performed using the full perturbation equation Eq. (f26|), including also 
the loop corrections. 



4 Results 

In this section we highlight the modifications to the tensor perturbation equations induced by 
the corrections adopted to regularise the background evolution. To this purpose, we shall first 
restrict the discussion to the tree-level a 1 corrections; the reason for this choice is twofold. 

First, a' corrections are derived unambiguously since we require them to reproduce the 
string scattering amplitude, as opposed to our lack of precise knowledge for quantum loop 
corrections, which forces us to speculate on plausible terms. Second, if any clear imprint of 
a "string phase" - like the one introduced in Jl7|, f|2|, |43| - is to be found in the spectrum of 
relic gravitons, it should be a direct consequence of a long enough string phaseQ characterised 
by a constant Hubble parameter. As suggested in []32] , this can be achieved by assuming that 
the evolution of the universe is driven by a 1 corrections into a fixed point in the (4>, H) space. 

In the final part of this section we will comment on the impact of the quantum loop 
corrections, and confront our predictions with the results obtained through exact numerical 
integrations of the perturbation equation. 



4.1 Tree-level a' corrections 

The lowest-order effective action of string theory provides an adequate description of the 
dilaton-driven branch of the pre-Big Bang scenario. In this weakly coupled and low curvature 
regime, no correction is required either to the effective potential nor to the wavelength of a 
given perturbation, and the evolution equation of each Fourier mode reduces to its usual form 



Eq. (|27| ). However, it is well known that the kinetic energy of the dilaton necessarily drives an 
isotropic and homogeneous background towards a curvature singularity, imposing the need to 
supplement the tree-level action with corrections in order to establish a reliable description of 
the high curvature regime. The very presence of the tree- level a' corrections Eq. (||), subject 
to the constraint Eq. (]?]), may lead in general to a regime of constant Hubble parameter and 
a linearly growing dilaton, H = const, <fi > 0, which seems to represent an exact conformal 



solution of string theory also to all orders in a' [[S^ • 

However, at any given finite order, we recall that not all sets of coefficients of the truncated 
action lead to such a good fixed point in the (0, H) plane (see for a detailed analysis). We 
shall thus restrict the analysis to the region of parameters which in principle may allow for 
a graceful exit, when loop corrections are included. In the (c 2 ,c 4 ) parameter space, such a 
region is delimited on the one hand by c 2 < 38/9, referring to a change of sign of the time 

derivative of the shifted dilaton, = <fi — 3H [f|7]]. On the other hand, tree-level a' corrections, 
truncated to first order, are not expected to violate the null energy condition, and thus cannot 



1 By long string phase, we mean N 3> 1, where N = ln(a) is the number of e-folds. 
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lead in principle to fixed points located after the Einstein bounce [f|7j (i.e, the bounce of the 
scale factor in the Einstein frame). This implies C4 > C2/2 for —00 < C2 < 4, to a first 
approximation. The location of such fixed points, satisfying the above constraints, can be 
seen in Fig. [I]. 




Figure 1: The figure on the left, from Ref. ||, shows the location in the {<p,H) space of 
the fixed points for different choices of the parameters C2 and C4, and the corresponding ratio 
£ = 4>/H. The blue line shows the C2-dependence of the fixed points for a constant C4 = 1, 
while the green line shows the C4-dependence for a constant C2 = 0. The point labeled by (1) 
corresponds to c 2 = 0, c 4 = 1, with ~ —0.72. The point labeled by (2) corresponds to 
c 2 = —40, c 4 = 1, with £( 2 ) — —0.56. The point labeled by (3) corresponds to c 2 = 0, c 4 = 15, 
with £( 3 ) ~ —0.13. On the right, we show the tree-level a 1 asymptotic shift c given by Eq. ( |35"|) 
as a function of the parameter C2 and of the fixed point distribution. The blue line shows the 
evolution of the shift for a constant C4 = 1. 



4.2 Effective frequency shift 

To describe the evolution of the high frequency modes, which are expected to leave the Hubble 
radius during the high curvature regime, we must use the perturbation equation Eq. (|26|), 
including those corrections required to regularise the background evolution. As mentioned 
in Section 3.2, the incorporation of the higher-order corrections in the perturbation equation 
implies a time-dependence of the effective wavenumber k. Indeed, Eq. (|26"|) can be conveniently 
rewritten 

M + {k 2 [l + c(ji)]-V(Ti)}l> h = Q, (34) 

where we have introduced the effective shift in the comoving frequency c = (y 2 — z 2 )/z 2 which, 
at next-to-leading order in the string tension expansion, reads 

_ 2a'{(4-c 2 )0' 2 -4(0"-2^)} 

c — -, — s . l^^J 

8a 2 - a'(j}' {8^ - c 2 (j)'\ 

In the weak coupling regime, and for small curvature (a' —>■ 0), the shift in the frequency 
is zero, hence no modification to the DDI cubic branch of the relic gravitational spectrum 
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is expected. However, this frequency shift will be no longer negligible when H ~ 0(a'). 
In general, we find that c(r]) tends to grow, although a short decrease (which still satisfies 
c{rj) > —1) may happen at the onset of the string phase. We also observe that the modification 
arising in the asymptotic regime, where the curvature is saturated by a linearly growing 



dilaton, corresponds to a constant shift of the comoving frequency, as first suggested in [[33 
In that case, by setting <fi = <f)'/a = <p and H = a' /a 2 = H, we can rewrite the asymptotic 
version of Eq. (^5]) as 

2a'0{(4-c 2 )0 + 4if} 



8 + a'4>{c 2 4>-8H} 



(36) 



During such an intermediate string phase the external potential reduces to its tree-level ver- 
sion, namely V(rf) = £ = a &~^^ 2 , and the enhancement in frequency associated to 
Eq. (^) is the same for all modes, since the shift c does not depend on the wavenumber k. 

In Fig. m we have represented the effective shift as a function of the parameter c 2 , high- 
lighting the C2-dependence by drawing the curve for C4 = 1. Although the asymptotic shift 
cannot be expressed as a function of an unique variable, it is however possible to extract its 
minimal value for the type of a 1 correction we are considering. The numerical simulations 
suggest that the lowest shifts correspond to fixed points located near the branch change line, 
(j) = <p — 3H = 0. The bound c 2 = 38/9 attached to the sign change of implies that, in 
the limit of large values of <j>, the effective asymptotic shift is always positive, and is given by 
c > 10/7. 

Does the frequency shift have any implication on the spectral distribution? In principle the 
answer is yes, since higher frequency modes tend to leave the Hubble radius at larger values of 
H in the pre-Big Bang scenario. As a consequence, a given comoving mode will spend less time 
under the potential barrier, resulting in a smaller amplification of the frequency attached to 
this mode. Although it is difficult to quantify such an effect, we believe the higher-curvature 
corrections amount to an overall rescaling, by a numerical factor of order unity, of the total 



energy density of the background, as already discussed in [33 



Finally, the form of the effective shift Eq. (|35|) suggests that a background undergoing 
a non-singular evolution in the vicinity of 8 — 8(f)H + c 2 2 = will automatically imply a 
rapid growth (if not a divergence) of the amplitude of the tensor fluctuations. Hence, we can 
understood the denominator of Eq. (|35| ) as a source of quantum gravitational instability, of 
the type already discussed in |48L [49[] . 



4.3 Full corrections and non-singular evolution 

In recent studies of the graceful exit problem, in the context of the pre-Big Bang scenario, the 
regularisation of the background curvature and the transition to the post-Big Bang regime, 
are obtained by including quantum loop corrections, according to Eq. (||). 

In this context, we may note that the tree- level a ; shift Eq. ( |35D provides an interesting 
constraint on the choice of parameters we can use in order to implement such a graceful exit, 
since in general the divergence curve for the high curvature shift is very close to the condition 
for the fixed point in the (0, H) plane. In 0, it has been shown that for a constant C4 the fixed 
point value of the Hubble parameter tends to be reduced when c 2 is decreasing. Simultaneously 
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the location of the singularity curve also decreases in the (0, H) plane, so that the range of 
coefficients leading to a non-singular evolution of the background, and avoiding the instability 
of quantum fluctuations, is strongly restricted to values very close to the simplest case of a' 
correction, C a i ~ {R%b ~ (^0) 4 } ( as a l so discussed in []32| ). This is an unexpected, although 
fortunate constraint on the loop corrections, which allows restricting the range of values for 
the parameters di and e;. 

In our attempt to characterise the power spectrum of metric perturbations generated in 
such a class of non-singular models, we will thus consider the perturbation equation Eq. ( p4|) 
with a full shift determined by 



2a' {(j)' 2 [4 - c 2 - Ad 2 e (f) - (4e x + e 2 )£e 2 ^ - 4 (j)" - 2^) [l + 4 J Be 1 e 2 *] } 



8a 2 - a'cj)' {8^ [1 + Be^} - 0' [c 2 + Ad 2 e^> + Be 2 e 2 <t>}} 



(37) 



Here the coefficients of the high curvature, one- and two loop corrections are respectively q, di 
and e,, with c\ = — 1 from the constraint Eq. (|7]). 

As already stressed, this frequency shift reduces to a constant in the asymptotic stage of the 
string phase. The loop corrections re-establish the time-dependence during the graceful exit, 
but its contribution soon becomes negligible at the end of the phase of accelerated evolution. 
Indeed, Eq. (|23|) - Eq. ( p4|) rapidly reduce to their tree- level expressions y ~ z ~ ae - ^ 2 , hence 
we expect k e ff — > k before we enter the radiation-dominated epoch, where the frequency shift 
vanishes identically. 

4.4 Observable implications 

We are interested in the generic features of the primordial gravitational waves, generated by 
quantum fluctuations of the background metric during the pre-Big Bang phase, which could be 
detected by both ground-based detectors such as LIGO |5D[ and VIRGO |5"I| , and space-based 



experiments such as LISA |52fl (see also |53fl for a recent review and references therein). For 
such a purpose, we use the dimensionless spectral energy density of perturbations, Q gw (uj), 
to describe the background of gravitational radiation. In critical units, the energy density is 
defined by 

^(w,^-^ (38) 
p c d In uj 

where dp gw is the present energy density in the stochastic gravitational waves per logarithmic 
interval of frequency, u(t) = k/a(t) is the physical (angular) frequency of the wave, and 
p c (t) = 3MpH 2 (t)/8ir is the critical energy density required to close the universe. 

Following the standard procedure for the computation of the energy spectrum ||45|| , we 
must determine the expectation number of gravitons per cell of the phase space, n(x, k) — rif, 
which only depends on the frequency for an isotropic and stochastic background. Neglecting 
a' and loop corrections, both in the background and perturbation equations, the lowest order 
result Eq. fl33|) leads to 

^(-,t)-^^(^) 3 , -<-i (39) 
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where u)\ = ki/a = H\a\ja is the maximal amplified proper frequency, associated with the 
value Hi of the Hubble expansion parameter at the end of the string phase. 

Including higher-order corrections, the relation between the ultraviolet cutoff of the spec- 
trum and the time of horizon crossing is in general more complicated, and a precise definition 
of the maximal amplified frequency is out of reach without a complete analytic solution for 
the background dynamics. However, a reliable estimate of the slope of the intermediate string 
branch of the spectrum can be obtained by using the leading features of the string phase. 

First, we shall neglect (during the string phase) the effective frequency shift in the pertur- 
bation equation, whose role is simply expected to induce an overall rescaling of the amplitude 



33fl , as already pointed out in Section 4.2. Second, we will assume that the duration of the 
"exit phase" is small compared to the intermediate string phase, and that no dramatic physics 
take place during the transition to the standard radiation-dominated epoch (indeed, in some 
cases where the exit is catalysed by the loop corrections, the process is almost instantaneous). 
This is a big assumption, of course, which could possibly underestimate "pre-heating" and "re- 
heating" effects. However, the accuracy of these assumptions can be confirmed with numerical 
simulations. 

In summary, we suppose that the evolution of the background consists mainly of three 
phases: an initial, low energy dilaton-driven branch, the intermediate string phase with a 
constant Hubble parameter (H) and a linearly growing dilaton (<f) = <fi = const), and the final 
radiation phase with decreasing Hubble parameter and frozen dilaton. It follows from our 
assumptions that the spectral graviton distribution can be correctly estimated through the 
low-energy perturbation equation, even for the high frequency branch of the spectrum. 



In that case, it is already well known [17] that the slope of the high-frequency modes, 



crossing the horizon in the high curvature, stringy regime, and re-entering in the radiation era, 
is fully determined by the fixed point values 0, H. Indeed, during such a string phase, the scale 
factor undergoes the usual de Sitter exponential expansion, while the logarithmic evolution of 
the dilaton, in conformal time, is weighted by the ratio <f)/H, i.e. <p(i]) ~ —(<p/H) log(— 77) + 

const fl54| . By introducing the convenient shifted variable (f> = <p — 3H, and referring the 

spectrum to a fixed point allowing a subsequent (loop catalysed) exit, i.e. < 0, one easily 
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fiflia ~ — , u 8 < uj < u x (40) 



\UJi. 

where uo s is the limit frequency marking the transition to the high curvature regime. 

In our model of the background evolution, on the other hand, the allowed fixed points <fr, H 
are located in the (4>,H) plane between the "Branch change line", <fi < 0, and the "Einstein 

Bounce", <fi < —H, see Fig. 1 and 0. It follows that, in the context of our approximations, 
the spectrum of tensor fluctuations, at high frequency, has a slope constrained by 

(Einstein Bounce) 2 < 3 — 4>/H < 3 (Branch change). (41) 

In order to check this important analytical result, we analysed numerical solutions for different 
choices of coefficients of the tree-level a' corrections. 
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4.5 Numerical results 



The spectral distribution of relic gravitational waves can be obtained by numerical integration 
of the perturbation equation Eq. (p6[) , in the background Eq. The initial conditions are 

chosen in the low curvature and weakly coupled regime and are thus based on the tree-level 
solutions, which provide an adequate description of this phase. We then evolve the system 
until a given mode re-enters the Hubble radius, in the late time FRW radiation-dominated 
epoch, where we extract the Bogoliubov coefficients by comparing the result of the simulation 
and the free-oscillating solution Eq. (fffi). 

To highlight the impact of the time- dependence of the frequency shift, we first present 
the spectral distribution in a case where we do not consider the a' and loop corrections in 
the perturbation equation Eq. ( j2"6|) . Fig. |2| illustrates the results of such a simulation,for the 
particular model of background corresponding to ci = d\ — e\ — — 1, C4 = d± — e± — 1 
(all the other d^d^a = 0), A = 1, B = -2 x KT 3 and = 5.63 x 10 -4 . We show, in 
particular, the non-singular evolution of H and <f), and the evolution of the quantity aH, 
which enables us to determine if a given mode leaves the Hubble radius during the DDI phase 
or during the intermediate string phase. Finally, we present the graviton distribution for such 
a background, computed from Eq. (|26|), in units k/k max , where k max = Max(aif) is a maximal 
amplified frequency. 




Figure 2: The choice of coefficients for the tree-level a' corrections is c\ — — 1, C2 = C3 = 
and C4 = 1 (and the same for di and e^), corresponding to the case (1) in Fig. |l|. The left figure 
shows a non-singular evolution for the Hubble parameter H = a/a and for 0/3, as a function 
of the number of e-folds, N = In a. The middle figure shows the evolution of aH /Max(aif) as 
function of N. The low-energy, dilaton-driven phase takes place approximatively for — 00 < 
N < —3. After a short transition, this initial period is followed by a string phase with 
nearly constant Hubble parameter and linearly growing dilaton, for 2 < iV < 9. After 
a successful exit triggered by loop corrections, the background evolution enters the FRW 
radiation-dominated phase at N ~ 16. Comoving modes leaving the Hubble radius during 
the string phase (10~ 6 < k/Ma,x(aH) < 1CT 2 ) are characterised by a slope smaller than 3 in 
the spectral distribution, as illustrated in the right figure which shows dpk/dhi(k) in units of 
k I k max - The upper dashed line corresponds to the cubic slope typical of the low frequency 
part of the spectrum, emerging from the dilaton-driven epoch. 

The effects of the higher-order corrections are expected to arise when the curvature scale 
becomes large in string units. As a consequence, the low frequency branch of the spectrum is 
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unaffected by such corrections, and remains characterised by a cubic slope. However, modes 
leaving the Hubble radius during the intermediate string phase are strongly affected by the 
background dynamics. The resulting slope of the spectrum is found to be reduced, but remains 
confined between 2 and 3, as expected from Eq. pip. 

Does the inclusion of the time-dependence of the comoving wavenumber lead to any modifi- 
cation? We stress that when such a shift in the frequency is considered, we have ky/1 + c = aH 
at horizon crossing. As a consequence, the shortest scale to leave the Hubble radius is given 
by k max = Max (aH/ \/T+~cV instead of M&x(aH). This is illustrated in Fig. H, where we 
compare the results of the numerical integration with and without corrections in the pertur- 
bation equation, for identical non-singular background evolutions (the same model of Fig. 
2). No modification appears during the low energy epoch, where the frequency shift is neg- 
ligible. However, higher-order corrections in the perturbation equation, which result in the 
time-dependent shift, lead to a reduction of the maximal amplified frequency and of the time 
spent on super-horizon scales. 




Figure 3: Using the same choice of coefficients as in Fig. |2|, we compare the evolution of 
aH /M&x(aH) (solid curve) and aH/Max (aH/ y/1 + cj (dashed curve) as a function of N. 
When higher-order corrections are included in the perturbation equation, the high frequency 
cutoff and the duration of the super-horizon epoch are found to be reduced. The right hand 
figure is a magnified image of the end point of the string branch of the spectral distribution, 
which highlights the impact due to a non-trivial frequency shift: the significant new feature 
is the reduction of the high frequency peak and a further reduction of the slope, which may 
also become smaller than two, for comoving modes leaving the Hubble radius during the "exit 
phase" . 

The impact of the time-dependence of the comoving wavenumber k is non-trivial for very 
high-frequency modes, leaving the Hubble radius during the graceful exit to the FRW-radiation 
dominated epoch. In that case, the result Eq.(|40D no longer applies, and the slope can be 
smaller than two, as shown in Fig. 3. However, the impact remains negligible for modes 
well inside the string branch of the spectrum i.e. for modes crossing the horizon during the 
asymptotic, fixed point regime. In that case, our analytical approximations for the slope of 
the spectrum are confirmed by the numerical simulations. In Table [lj we present the predicted 
and measured slopes well inside the string phase, for the three different choices of coefficients 
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illustrated in Fig. [l|. The agreement is striking. 





Coefficients 


at O(a') 


Predicted slope 


Measured slope 
(numerical) 


c 2 


c 4 


3- 


4>/H 




Case (1) 







1 


2.44 


2.45 


Case (2) 


-40 




1 


2.28 


2.3 


Case (3) 







15 


2 


87 


2.85 



Table 1: Choices of coefficients for the tree-level a 1 corrections with c% = — 1 and the remaining 
coefficient satisfying c 3 = — (ci + c 4 + c 2 /2), according to Eq. (0). We have compared the 
predicted slope determined by the fixed point, according to Eq. (ffij), with the measured slope 
well inside the string branch of the spectrum, obtained by numerical integration of Eq. fl26|) . 
including all corrections. 



It is well known |44[] that the spectrum of relic gravitons, in the context of the pre-Big 
Bang scenario, easily evades constraints arising from both the CMB anisotropy at COBE scale 



55] and pulsar timing data [pq| , because of its rapid growth at low frequency. Hence, the 
peak amplitude of the spectrum is only constrained by the primordial nucleosynthesis bound 
|57|, as well as by primordial black hole production |58| . 

The relic graviton background is in general compatible with these constraints if it is nor- 
malised so as to match the string scale at the high-frequency end point of the spectrum. In 
that case one finds j^TJ that the high frequency peak is typically of order Vt gw {oj) ~ 1CT 6 , 
for a maximal amplified frequency u ~ 10 11 [Hz]. Saturating this high-frequency end point, 
and using the prediction Eq. ([IT]) for the slope of the string branch of the spectrum, the 
energy density of relic gravitons from a pre-Big Bang phase could be at most of O(l0~ 15 ) at 
uj ~ 10 2 [Hz], regardless of the duration of the string phase with constant Hubble parameter, 
i.e. far below the sensitivity of the second (planned) generation of interferometric gravity wave 
detectors . It seems thus difficult to detect the relic gravitons from a high curvature string 
phase, associated to a fixed point of the truncated effective action - at least within the model 
of background considered in this paper - unless the exit phase is able to affect a sufficiently 
wide band of the high frequency spectrum. 

Finally, the numerical integration of Eq. (^) enables us to follow carefully the time evo- 
lution of modes, from the DDI epoch up to the FRW-radiation dominated era. We have 
explicitly checked that the amplitude of tensor perturbations is rapidly frozen in on super- 
horizon scale, regardless of the particular higher-order dynamics of the background during the 
string phase and the exit epoch. 



5 Conclusion 

In this paper we have discussed the evolution of tensor perturbations in a class of non-singular 
cosmological models based on the gravi-dilaton string effective action, expanded up to first 
order in a', and including one-loop and two-loop corrections in the string coupling parameter. 
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The coefficients of the loop expansion have been chosen in such a way as to avoid the 
curvature singularities in the background, and the quantum instability of metric fluctuations. 
In addition, we have included an appropriate radiation backreaction to stabilise the dilaton in 
the final, post-Big Bang regime. We have been able to obtain, in this way, a class of regular 
cosmological solutions in which the initial Universe, emerging from the string perturbative 
vacuum, is first trapped in a fixed point of high, constant curvature, and then spontaneously 
undergoes a smooth transition to the standard FRW regime. The fixed point is determined 
by the a' corrections, while the unstability and the decay of the high curvature regime is 
triggered by the quantum loop corrections. 

We have determined, in this class of backgrounds, the linearised equation governing the 
evolution of tensor metric fluctuations, taking into account all the a' and loop corrections 
contributing to the background regularisation. We have discussed, on this ground, the ampli- 
fication of tensor perturbations normalised to the quantum fluctuations of the vacuum, and 
we have estimated their final spectral energy distribution through analytical and numerical 
methods. 

Our results confirm previous expectations, that the low frequency modes, crossing the 
horizon in the low-curvature regime, are unaffected by higher-order corrections. Also, we have 
explicitly checked (through numerical simulations) that the so-called "string branch" of the 
spectrum, associated with the asymptotic, fixed point regime, can be reliably estimated by 
means of the low-energy perturbation equation, that its slope is determined by the asymptotic 
constant values of H and <fi, and that it is flatter than at low frequency. 

However, when the background regularisation crosses a fixed point, and when the fixed 
point is determined by the a' expansion truncated to first order (like in the class of models 
of this paper), it turns out that the slope of the string branch is at least quadratic, and thus 
probably too steep to be compatible with future observations by planned advanced detectors. 
In the transition regime dominated by the quantum loops, on the other hand, the influence of 
the higher-order corrections on tensor perturbations is more dramatic, and the corresponding 
slope may be much flatter than the one of the preceding string phase. 

In summary, the results of this paper agree with the general expectation that the shape 
of the spectrum of the relic graviton background, obtained in the context of the pre-Big 
Bang scenario, is strongly model-dependent. However, it turns out that the formulation of a 
complete and consistent regular scenario, associated to a "visible" graviton spectrum, may be 
more difficult than previously expected - at least within a truncated perturbative approach, 
like the one adopted in this paper. 
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